LA-UR-02-7498 



Status of Bk from Lattice QCD 



Raj an Gupta 1 % 

Theoretical Division, Los Alamos National Lab, 
Los Alamos, New Mexico 87545, USA 



Abstract 

A brief review of lattice calculations of the bag parameter Bk relevant for understanding indirect 
CP violation in the neutral kaon sector is given. A status report on current state-of-the-art calcula- 
tions is presented as well as a discussion of the value of Bk exported to phenomenologists. This is a 
condensed and updated version of the review presented at the CKM Unitarity Triangle Workshop 
held at CERN during Feburary 13-16, 2002 (http://ckm-workshop.web.cern.ch/ckm-workshop/). 
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I. INTRODUCTION 



The most commonly used method to calculate the matrix element 

K° | Z (sd)v-A(sd)v-A(n) | K ^ is to evaluate the three point correlation function 

shown in Fig. ^ This corresponds to creating a K° at some time t\ using a zero-momentum 
source; allowing it to propagate for time to — ti to isolate the lowest state; inserting the 
four-fermion operator at time to to convert the K° to a K°; and finally allowing the K° to 
propagate for long time t 2 — t - To cancel the K° (K°) source normalization at times t\ 
and t 2 and the time evolution factors e~ ERt for times t 2 — to and to — ti it is customary to 
divide this three-point function by the product of two 2-point functions as shown in Fig 1. 
If, in the 2-point functions, the bilinear operator used to annihilate (create) the K° (K°) 
at time to is the axial density S74750?, then the ratio of the 3-point correlation function to 
the two 2-point functions is (S/3)Bk- 

Bk is defined to be the value of the matrix element at the physical kaon and normalized 
by the Vacuum Saturation Approximation value 8/3M K F K 

(K G I Z {Hd) v - A {sd)v-A{n) I K 5 ) = (S/3)B K M 2 K F 2 K . 

Earliest calculations of Bk were done using Wilson fermions and showed significant devia- 
tions from this behavior. It was soon recognized that these lattice artifacts are due to the 
explicit breaking of chiral symmetry in the Wilson formulation 0, 0, H 0, |H[ . Until 1998, the 
only formulation that preserved sufficient chiral symmetry to give the right chiral behavior 
was Staggered fermions. First calculations using this approach in 1989 gave the quenched 
estimate Bx(NDR,2GeV) = 0.70 ± 0.01 ± 0.03. In hindsight, the error estimates were 
highly optimistic, however, the central value was only 10% off the current best estimate, 
and most of this difference was due to the unresolved 0(a 2 ) discretization errors. 

In 1997, the staggered collaboration refined its calculation and obtained 0.62(2) (2) 
again the error estimate was optimistic as a number of systematic effects were not fully 
included. The state-of-the-art quenched calculation using Staggered fermions was done by 
the JLQCD collaboration in 1997 and gave B K (2GeV) = 0.63 ± 0.04 @. This estimate 
was obtained using six values of the lattice spacing between 0.15 — 0.04 fermi, thus allowing 
much better control over the continuum extrapolation as shown in Fig. |2] along with other 
published results. This is still the benchmark against which all results are evaluated and is 
the value exported to phenomenologists. This result has three limitations: (i) It is in the 
quenched approximation, (ii) All quenched calculations use kaons composed of two quarks 
of roughly half the "strange" quark mass and the final value is obtained by interpolation to a 
kaon made up of (m s /2, m s /2) instead of the physical point (m s , m^). Thus, SU(3) breaking 
effects (m s 7^ m^) have not been incorporated, (iii) There are large 0(a 2 ) discretization 
artifacts, both for a given transcription of the AS* = 2 operator on the lattice and for different 
transcriptions at a given value of the lattice spacing, so extrapolation to the continuum limit 
is not as robust as one would like. These limitations are discussed after a brief summary of 
the recent work. 

In the last four years a number of new methods have been developed and the correspond- 
ing results are summarized in Table 1. 

• The Rome collaboration has shown that the correct chiral behavior can be ob- 
tained using 0(a) improved Wilson fermions provided non-perturbative renormal- 
ization constants are used. Their latest results, with two different "operators", are 
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FIG. 1: 



Collaboration 


year 


B K (2GeV) 


Formulation 


Renormalization 


a" 1 (GeV) 


Staggered [6] 


1997 


0.62(2)(2) 


staggered 


1-loop 


oo 


JLQCD [7J 


1997 


0.63(4) 


staggered 


1-loop 


oo 


Rome [8] 


2002 


0.63(10) 


Improved Wilson 


NP 


oo 


Rome [8] 


2002 


0.70(12) 


Improved Wilson 


NP 


oo 


CP-PACS [9] 


2001 


0.58(1) 


Domain Wall 


1-loop 


1.8 GeV 


CP-PACS [9J 


2001 


0.57(1) 


Domain Wall 


1-loop 


2.8 GeV 


RBC [10] 


2002 


0.53(1) 


Domain Wall 


NP 


1.9 GeV 


DeGrand [11] 


2002 


0.66(3) 


Overlap 


1-loop 


1.6 GeV 


DeGrand [11] 


2002 


0.66(4) 


Overlap 


1-loop 


2.2 GeV 


GGHLR [12] 


2002 


0.61(7) 


Overlap 


NP 


2.1 GeV 



TABLE I: Quenched estimates for Bk evaluated in the NDR scheme at 2GeV. The fermion 
formulation used in the calculation, the method used for renormalizing the operators, and the 
lattice scale at which the calculation was done are also given. NP indicates non-perturbative 
renormalization using the RI/MOM scheme and a -1 = oo implies that the quoted result is after a 
continuum extrapolation. 



B K (2GeV) = 0.63(10) and 0.70(12) |g. These, while demonstrating the efficacy of 
this method, do not supplant the staggered result, as the continuum extrapolation is 
based on only three points and the data have larger errors. The discretization errors 
can be characterized as Bx(a) = Bk(1 + a A) with A « 400 MeV and are similar in 
magnitude to those with staggered fermions at l/a — 2 GeV, as are the differences in 
estimates with using different operators. In the staggered formulation, the artifacts 
are, however, 0(a 2 A 2 ) and 0(a 2 ) and the data suggest an unexpectedly large A ~ 900 
Mev. 

Four collaborations have new results using domain wall and overlap fermions as shown 
in 

Table HI 0,03, El E3, 1 

Hill ]. Both formulations have built in chiral symmetry at finite 
a and 0(a) improvement. Each of these collaborations have used slightly different 
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FIG. 2: Published estimates of Bk with fermion formulations that respect chiral symmetry. All 
results are in the quenched approximation. 



methodology, so they cannot be compared head on, or combined to do a continuum 
extrapolation. Thus, the results are quoted with reference to the lattice spacing at 
which the calculation was done. The differences reflect 0(a 2 ) (and 0(a 2 ) in cases 
where perturbative renormalization constants have been used) artifacts. 

• Calculations using another method with good chiral behavior, twisted mass QCD, are 
in progress [HI]. 

Starting with the current best quenched lattice estimate, the JLQCD staggered result 
£>x(2GeV) = 0.63(4), deriving an estimate for the physical B K requires consideration of the 
following issues. 

• The 0(a 2 ) errors in the staggered formulation are large. Nevertheless, the error 0.04 
obtained by the JLQCD collaboration on including both 0(a 2 ) and 0(a 2 ) terms in the 
extrapolation is a reasonable la estimate of both the statistical and the extrapolation 
to continuum limit errors. 

• A choice for a s and the number of flavors in the perturbative expression has to be 
made to convert Bk — > Bk- It turns out that the result is insensitive to whether 
one uses quenched or full QCD values. Using the 2-loop expression, the result for the 
central value is B K = 0.86(6). 

• An estimate of the systematic uncertainty associated with the quenched approximation 
and SU(3) breaking. Preliminary numerical estimates suggest that dynamical quarks 
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would increase the value by about 5% [p| lil . Sharpe estimates, using ChPT, that 
unquenching would increase Bk by 1.05 ± 0.15, and SU(3) breaking effects would also 
increase it by 1.05 ± 0.05 ^t|- This analysis of systematic errors is not robust and, 
furthermore, the two uncertainties are not totally independent. So one can take an 
aggressive and a conservative approach when quoting the final result for Bk- In the 
aggressive approach, the error estimate is given by combining the central values in 
quadratures. This gives a 7% uncertainty and 

B K = 0.86 ± 0.06 ± 0.06 . (1) 

In the conservative approach, advocated by Sharpe jUJ , one combines the uncertainty 
in quadratures to get a 16% uncertainty. The final result in this case is 

B K = 0.86 ±0.06 ±0.14. (2) 

Given the lack of a robust determination of the systematic error, it is important to decide 
how to fold these errors in a phenomenological analysis. The recommendation is to assume a 
flat distribution for the systematic error and add to it a gaussian distribution with o = 0.06 
on either end, and do a separate analysis for the aggressive and conservative estimates. In 
other words, a flat distribution between 0.72 — 1.0 for a conservative estimate of Bk (or 
0.80 — 0.92 for the aggressive estimate) to account for systematic errors due to quenching 
and SU(3) breaking. Since this is the largest uncertainty, current calculations are focused 
on reducing it. 

Finally, the reasons why the quenched lattice estimate of Bk has been stable over time 
and considered reliable within the error estimates quoted above are worth reemphasizing: 

• The numerical signal is clean and accurate results are obtained with a statistical sample 
of even 50 decorrelated lattices. 

• Finite size effects for quark masses > m s /2 are insignificant compared to statistical 
errors once the quenched lattices are larger the 2 fermi. 

• In lattice formulations with chiral symmetry, the renormalization constant connecting 
the lattice and continuum schemes is small (< 15%), and reasonably well estimated 
by one-loop perturbation theory. 

• For degenerate quarks, the chiral expansion for the matrix element has no singular 
quenched logarithms (they cancel between the AA and VV terms) that produce large 
artifacts at small quark masses in observables like M%, f n , etc. Also, the chiral expan- 
sion between the quenched and full theories have the same form EaHHEl- 



ChPT estimates of quenching and SU(3) breaking systematic errors are at the 7—16% 
level [HElIIU. 
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